Dynamics of a colloid-stabilized cream 
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We use x-ray photon correlation spectroscopy to investigate the dynamics of a high volume fraction 
emulsion creaming under gravity. The dodecane-in-water emulsion has interfaces stabilized solely by 
colloidal particles (silica). The samples were observed soon after mixing: as the emulsion becomes 
compact we discern two regimes of ageing with a cross-over between them. The young emulsion 
has faster dynamics associated with creaming in a crowded environment accompanied by local 
rearrangements. The dynamics slow down for the older emulsion although our studies show that 
motion is associated with large intermittent events. The relaxation rate, as seen from the intensity 
autocorrelation function, depends linearly on the wave vector at all times; however, the exponent 
associated with the line shape changes from 1.5 for young samples to less than 1 as the emulsion ages. 
The combination of ballisticlike dynamics, an exponent that drops below 1 and large intermittent 
fluctuations has not been reported before. 

PACS numbers: 82.70.Dd, 82.70.Kj, 61.05.cf 



I. INTRODUCTION 

Colloid-stabilized emulsions are metastable composites 
which can be remarkably robust. The interfacial col- 
loids are trapped extremely strongly [l| because they re- 
duce the shared area between the dispersed and contin- 
uous phases. When an emulsion is initially prepared the 
droplets coalesce at a rate that decreases with increas- 
ing size corresponding to an increasing area fraction of 
interfacial colloids Ultimately a complete coverage 
of the interface by the particles provides a steric barrier 
against coalescence (although other mechanisms of sta- 
bilization are also possible Q [1] Hi)- Unlike the inter- 
faces of a conventional, surfactant stabilized, emulsions 
the densely packed colloids create a semi-solid layer 6]. 
This means that anisotropic stresses can be supported 
by these liquid-liquid interfaces 0] and that a high 
volume fraction emulsion is a novel solid Q . 

Here we report an experiment on colloid-stabilized 
emulsions as they become a compact cream layer; only 
preliminary studies currently exist fllj . In our sam- 
ples the typical final droplet sizes are on the scale of ten or 
more microns. The coalescence rates observed by Arditty 
and coworkers ^] strongly suggest that by the time our 
measurements begin the droplets have already reached 
full and stable coverage. Hence, this is a non-thermal 
system where gravity will drive robust droplets to form 
a high volume fraction cream: a process of solidification. 
For our preparation route there is a broad distribution of 
droplet sizes which implies that the cream will not tend 



to crystallize as the volume fraction increases: the likely 
end state of this system is a disordered random packing 
of droplets. The distribution also implies that there will 
be a spread of creaming rates that will necessarily lead 
to collisions between droplets from the earliest times. As 
described below, the emulsion begins at a volume fraction 
of = 0.4 and then becomes increasingly compact; from 
a very early stage during the experiment all the droplets 
are in contact. 

The creaming of conventional emulsions has been 
probed in great detail because it is of both fundamen- 
tal and industrial importance [12,] . In conventional emul- 
sions three generic types of behavior have been identified: 
I. If the emulsion is monodisperse with no inter-droplet 
attractions a close-packed cream layer can form while a 
clear solvent layer appears at the bottom. Between these 
two phases is a low volume fraction dispersion which is 
eradicated as the lower and upper interfaces move to- 
wards each other during creaming. II. If the emulsion 
is polydisperse the lower interface becomes diffuse rather 
than sharp while the volume fraction of the cream layer 
continuously increases. HI. If the droplets have attractive 
interactions a gel can form which may slowly compactify 
or suddenly collapse. No low volume fraction dispersed 
phase is observed. In all three cases the compact cream 
layer can have some similarities to a liquid-liquid foam. 
The coarsening of foams has been studied using diffusive 
wave spectroscopy and it has been observed that ripen- 
ing of bubbles leads to the build-upof stress followed by 
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sudden avalanche rearrangements 

Our samples are opaque to visible light and contain 
constituents with high electron density contrast. In ad- 
dition the dynamics are slow and we do not know prior to 
the experiment whether the samples will be ergodic. To 
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FIG. 1: (Color online) (a) The form factor of the particles used in this experiment measured using small-angle x-ray scattering; 
(b) raw x-ray data showing a snap-shot of a speckle pattern; (c) confocal microscopy image of a particle-stabilized emulsion. 
Dodecane is the dispersed phase; the white layers are the silica colloids trapped at the interfaces. The scale bar is 50 fim; (d) 
photographs of an emulsion ($ = 0.25 to ease visualization) illustrating creaming taking place. 



explore these systems a powerful approach is to use X-ray 
Photon Correlation Spectroscopy (XPCS) using a charge- 
coupled detection device (CCD) as a two-dimensional de- 
tector. The samples are in the single-scattering regime 
for 8 keV x-rays and, by recording the scattering at many 
pixels corresponding to the same wave vector, ensemble 
averaging can be carried out. Recently developed tech- 
niques can be used to reveal a detailed description of 
the fluctuations of the sample about its average behav- 
ior di 

Recent research on some soft-glassy materials using 
photon correlation spectroscopy both with visible light 
and with x-rays has shown a linear relationship between 
the probed length scale and the relaxation time in combi- 
nation with faster than exponential relaxation. Further- 
more, the dynamics are often characterized by intermit- 
tent events and show a delicate dependency on sample 
age This behavior is in stark contrast to Brown- 

ian dynamics but its origin is so far not fully understood. 
While a few soft glassy materials exhibit intermittent dy- 
namics as a result of small temperature fluctuations [l7[ 
or for undiagnosed reasons [l^ most systems are respond- 
ing to some driving force. An emulsion creaming under 
gravity is just such a system; again we find that the late- 
stage macroscopically quiescent system is remarkably ac- 
tive on a microscopic scale. 

II. MATERIALS AND METHODS 

The liquids forming the emulsion are water (Elga Pure- 
lab, 18.2 MOcm) and dodecane (Aldrich, >98%) used 
without further purification. The anticipated interfa- 
cial tension of 72.8 mN/m at 20° C ^ leads to ex- 
tremely strong trapping of the colloids used (more than 
8 X 10^ k^T per particle if the wetting angle of our 
particles lies between 60° to 120°). Water-oil interfaces 
can become charged during emulsion preparation 20] and 
this can inhibit the trapping of colloids on interfaces. To 
suppress this effect we add KCl (Aldrich, 99+%) suf- 
ficient to give a Debye screening length of 0.2 nm. The 
particles are silica prepared using the Stober method [2l[ . 
The particle radius, 243 ± 7 nm, was determined by fit- 



ting the form factor for a sphere to the time-averaged 
small-angle x-ray data [12] as shown in Fig.[lja). An ul- 
trasound probe (Sonics & Materials) is used to disperse 
the particles, 0.110 g (5vol%) or 0.165 g (7.5vol%), in 
0.6 ml water. Then the dodecane (0.4 ml) and KCl 
(0.1 mg) are added before further mixing. A syringe with 
long needle was used to fill a glass capillary (Hilgenberg) . 
Some of the same sample was also examined using bright- 
field microscopy to confirm that an emulsion formed. 

The experiment was carried out at the Troika beam- 
line at the ESRF with an 7.979 keV x-ray beam (A = 
0.155 nm). We used a channel cut Si (111) monochro- 
mator (AA/A « 10~^) and a 10 /im x 10 /im aperture 
before the sample was used to create a partially coherent 
beam. The scattering from the colloid-stabilized emul- 
sion was recorded using a CCD with 22.5 /im pixel size 
located 2.3m after the sample (see Fig. [TJb)) covering 
wave vectors 0.01 nm~^ < q < 0.08 nm~^. Analyz- 
ing the intensity variations yields quantitative informa- 
tion about the dynamics occurring within the sample. 
This is found from the intensity autocorrelation function: 
52(9, t) (/(g, 0)/(g, r))/(/(gf))^ where (...) denotes the 
average over a ^^20 minute time window and a ring of 
pixels sharing the same wave vector, q. The intermediate 
scattering function f{q,T) is proportional to the square 
root of 52(9, ■'■) ~ 1; we model this using a modified ex- 
ponential function: 

g2iq,T)-l = b + ce^p(^-2(^^y^ (1) 

with b as the value to which the curve decays; c is the 
apparent contrast which defines the difference in height 
between the initial and final level of the curve and Tc is 
the relaxation time ,16i] . The function 32(9, t) — 1 can 
be exponential (/3=1), compressed exponential (/3 >1) 
or stretched exponential (/3 <1). Extracting parameters 
giving the best-fit between the model and the data for 
different q allows us to establish a relationship between 
the length scale (i.e. inversely proportional to the wave 
vector) , the characteristic relaxation time Tj, and the ex- 
ponent p. 

We examine the evolution of the intensity autocorrela- 
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FIG. 2: (Color online) Intensity autocorrelation function calculated at q — 0.033 nm ^ for two different ages as written on the 
graphs: early (squares) and late (triangles). Colloid volume fractions (a) 5% (b) 7.5%. 



tion function over time to establish whether ageing and 
intermittency are present. To do so the degree of corre- 
lation, ci{t,T) is calculated for a set of fixed lag times r 
at a given wave vector q • 



Cl{t,T) 



{Ip{t)Ipit,T))p 
{Ip{t))p{Ij,{t,T))p 



(2) 



where Ip denotes the intensity at pixel p and r is the lag 
time between two frames while (. . .)p denotes the aver- 
age over all pixels contained in a given annulus of the 
CCD image corresponding to a given q. It is possible 
to detect intermittent behavior in the ci(t,T) signal by 
examining the variation of ci{t,T) around its mean for 
different lag times fl5|. Homogeneous dynamics show a 
variance of ci(t,T) that varies with ci{t^T); however, for 
heterogeneous dynamics the variance will show a peak 
around the typical relaxation time i.e. the time when 
most of the intermittent events take place. This method 
will only resolve such a peak around Tc if the fluctuations 
of the noise are smaller than the signal due to the particle 
dynamics. 



III. RESULTS AND ANALYSIS 

Here we compare the dynamics and ageing of two oil- 
in- water emulsions containing = 5% and 7.5% of silica 
colloids. As observed by confocal microscopy the typi- 
cal droplet sizes are 16 ± 4 fim and 12 ± 3 /im diame- 
ter respectively, e.g. Fig. [TJc). Measurements began as 
soon after emulsification as possible. The time-resolved 
speckle patterns were recorded for over 2 hours (with 
the samples returned later for subsequent measurements) 
and from these images the time-resolved correlation func- 
tion ci{t,T) is calculated, Eq. 



Figure [5] shows the change in g2{T) — 1 with age for the 
two different samples a.t q — 0.033 nm^^. A structure fac- 
tor peak due to static particle-particle correlations would 
occur a.t q < 0.013 nm"'^ which is at the lower limit of 
our wave vector range (marked by the dashed vertical 
line in Fig. [3]). Hence we are exploring the dynamics at 
or smaller than the particle scale. 

The emulsion behavior can be divided into two regimes 
with a cross-over between them. Initially the relaxation 
time is short, then there is a transition where the system 
starts slowing down and in the final regime it reaches a 
significantly longer relaxation time. For the 5% sample 
(larger droplets) the initial decay is significantly faster 
than that of the 7.5% sample (smaller droplets). The on- 
set of the slowing down occurs in both cases after about 
1 hour 20 minutes but is drawn out over almost an hour 
for the small droplets while the change from fast to slow 
relaxation happens over 10 minutes for the 5% sample. 
The final slow regime in the 5% sample reaches an ini- 
tial plateau; however, significant fluctuations in the re- 
laxation time are observed over the remaining time for 
which measurements were made. The slower ageing evo- 
lution of the 7.5% sample leads to only a short segment 
of the final stage being observed; however, this still seems 
to plateau at a similar relaxation time as the 5% sample. 

Shown in Fig.[ljd) is a more dilute emulsion ($ — 0.25) 
before and after creaming for one hour; the cream layer 
occupies well over half of the sample volume. For the 
x-ray studies, $ = 0.4, giving a cream layer that occu- 
pies most of the sample volume: the x-ray beam passes 
through the cream layer at all times and for all samples. 
Repeat measurements were made after the sample cap- 
illary had been removed from the instrument and later 
remounted (data not shown). In some cases the distur- 
bance to a late stage sample led to a return to early stage 
behavior even though the sample itself was then hours 
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FIG. 3: (Color online) (a, b) Relaxation rate 1/tc versus wave vector (c, d) exponent (3 versus wave vector, see Eq. ([T}: early 
(squares) and late (triangles). The dashed vertical line corresponds to the wave vector q = 0.013 nm~^ which corresponds to 
2n/d, where d is the diameter of the silica particles, (a, c) 5%, (b, d) 7.5%. Inset to (b): Relaxation rate at late times for both 
samples with an expanded y-scale. 



or days old. Since the coalescence of colloid-stabilized 
droplets is completed within minutes 0] we do not asso- 
ciate the early stage regime with coalescence. Instead it 
appears likely that early stage behavior can be a conse- 
quence of disrupting the droplets as an inevitable part of 
mounting the sample at the beamline. 

Qualitatively, the intensity autocorrelation function 
(Fig.mfor q = 0.033 nm~^) shows a significant increase in 
the relaxation time with age for both volume fractions. 
At the same time it is apparent that the shape of the 
decorrelation changes while the sample ages. To quan- 
tify these effects the best fit of Eq. ([T]) to all the obtained 
curves was found. The optimal fit parameters for all wave 
vectors are shown in Fig. [3) For diffusive dynamics we 
expect the relaxation rate 1/tc ~ and f3 — 1; for bal- 
listic motion we expect 1 /tc ~ q with a value of /3 which 
reflects the velocity distribution [l^l (/? = 2 corresponds 
to a Gaussian distribution). Evidently, 1/tc varies lin- 



early with q for both samples and at all times showing 
that the dynamics leading to the relaxation are not dif- 
fusive; however, there is a significant change between the 
behavior at early and late times. For early times there is 
a sizeable off-set such that 1 /tc does not decrease to zero 
for small q. This is observed for both volume fractions 
and is also evident in other data sets not presented here. 
In addition to this off-set the slope of the 1 / Tc-curve de- 
creases with time as the droplets are compactified (the 
slope is smaller for the smaller droplets, Fig.[3Ib) inset). 

The background level, 6, to which the intensity auto- 
correlation function decays, eq. ([1]), shows slight wave- 
vector dependence. We observe non-zero values for 
.92(9, r) — 1 at long times for q values close to the maxima 
of the form factor of the particles. Non-zero g2{q, ''') — 1 
values at long times are typically associated with non- 
ergodic sample dynamics; however, associating error bars 
with these long time measurements is challenging. In the 
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FIG. 4: (Color online) The signal ci{t,T) at q = 0.033 ± 0.003 nm ^ for different r values plotted against sample age shows 
that large fluctuations are present, (a) 5%, (b) 7.5%, lag times given for curves from top to bottom. 



early stage the existence of a component of 1/tc which 
is independent of wave-vector impUes that the samples 
have dynamics that lead to complete decorrelation at a 
steady rate. The long-time behavior is hence not due to 
non-ergodicity but instead appears to be a result of the 
small size of the signal. For late stage creaming behavior, 
where the dynamics are slow, there are insufficient data 
to firmly determine the long-time behavior of the sample. 

Although the time scales change for emulsions stabi- 
lized using two different volume fractions of colloids the 
overall aging behavior does not. At early times the ex- 
ponent (3 has values around 1.5 for smaller wave vec- 
tors probed (see Fig. ^c) and (d)) and slowly decreases 
towards 1.0 for the larger wave vectors [l^. At later 
times P is smaller starting at 1.0 for lower wave vec- 
tors and decreasing further for the largest wave vec- 
tors probed. The values imply that relaxation is faster 
than exponential (compressed) at least at early times. 
Such behavior, which has also been found by numerous 
other researchers in a wide variety of samples, cannot 
be created from the sum of exponential processes and 
its interpretation is the subject_of great current activ- 
ity (a 
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times the intensity autocorrelation function is a slightly 
stretched exponential: again this is not typically seen in 
concert with ballisticlike relaxation rate behavior. 

To determine whether the average dynamics capture 
what is essential about this system we need to establish 
if the rearrangements under gravity are continuous over 
time, representing a slow pushing of droplets together, 
or instead the upward drive induces a multitude of sepa- 
rate, small rearrangements. To carry out this analysis an 
optimized g-ring (large number of pixels, high intensity) 
was chosen and time sequences are selected (as described 
below) to calculate the variance of ci{t, r) over a specific 
time range for the measured r-range. Figure 2] shows the 
fluctuations of ci{t,T) for different values of r over the 
measurement time of the sample. The two regimes of 



sample dynamics are clear in these raw data. B y c alcu- 
lating cj(t, 0) which by definition is a static value |l5| and 
examining its fluctuations it is evident that noise forms a 
substantial fraction of some time sequences. It is not pos- 
sible to directly correct for this noise using the number of 
pixels available. Instead, segments of ^^20 minute dura- 
tion, in which the noise seen in c/(i,0) is less significant, 
are chosen for determining the variance of the c/(i,T)- 
signal. A large variance will indicate that an important 
role is being played by more extreme events. 

Figure [5{a, b) shows the time averaged correlation 
functions g2{q,T) — 1 for 5% and 7.5% with the vari- 
ance of the underlying c/(t, r) signal over time t at each 
T overlaid. Especially for the late 5% sample a clear peak 
around the relaxation time is visible suggesting that in- 
termittent rearrangements occur at late stages. Less of 
the late stage was accessible for the 7.5% sample. For 
the early times (data not shown) it does not appear that 
the variance shows a peak but due to the lack of data 
for small enough r this cannot be said for certain. There 
is no indication that these intermittent events are asso- 
ciated with a collapse of the emulsion; it is most likely 
that our emulsions are type II. 



IV. DISCUSSION 

In these experiments the x-rays are scattered from 
colloidal silica, irreversibly trapped at the surface of 
oil droplets. The experimentally accessible wave vector 
range corresponds to the colloidal scale rather than the 
droplet scale. Hence our samples are organized into large 
units with correlated dynamics on the scale we are ob- 
serving. The droplets are driven to rise and pack more 
densely by gravity and the colloids are carried along as 
this process takes place. 

Because a linear variation of relaxation rate with wave 
vector is a signature of ballisticlike dynamics we need to 
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FIG. 5: (Color online) A peak in the variance in the region of the relaxation time is an indication for intermittent behavior. 
Triangles: c/(t,r) for optimized g-range and time segment (late stage) that does not show strong photon fluctuations, circles: 
variance of c/(t,r) for different r. (a) 5% and (b) 7.5%. 



begin by considering how fast we would expect creaming 
to occur. Yan and Masliyah have applied the concepts of 
hindered creaming to particle-stabilized emulsions p^ . 
The analysis is based on the model of hindered settling 
due to Richardson and Zaki: the creaming velocity is 
expected to be Wc = wo(l — ^)"' where vq is the creaming 
velocity of the unhindered subunits and the exponent n 
depends on the interactions between the subunits (with a 
value around 5 being typical) ^ [30,]. For our droplets 
at volume fraction (f> = 0.4, using n = 5, this analysis 
gives a hindered creaming velocity of Vc ~ 300 nm/s. 

The early stage creaming under gravity is character- 
ized by a relaxation that is independent of wave vector 
as indicated by the intersection of the 1 / Tc-line with the 
y-axis (see Fig. ^a) and (b)). From the intercept we 
extract the rate of 'whole sample changing' events. The 
5% and 7.5% samples fully decorrelate at quite a high 
rate (every 14 and 60 seconds respectively at this age). 
The time it takes for the creaming sample to pass through 
the 10 ^m beam (about 30 sec) according to the Richard- 
son and Zaki model is roughly right to account for these 
decorrelations. The difference in the rate of 'whole sam- 
ple changing' events between the two emulsion composi- 
tions is larger than predicted by the model of hindered 
settling. We note that the hindered creaming velocity 
depends on the volume fraction of the emulsion to a high 
power (here taken to be 5). Small changes in the compo- 
sition of the emulsion transferred to the x-ray capillary 
will be reflected in large changes in the creaming velocity. 
It is likely that this accounts for the differences observed. 

Upward sample motion accounts for the wave vector 
independent aspect of the relaxation rate; the linear q- 
dependent relaxation process rides on the back of the 
steady creaming and is presumably associated with lo- 
cal rearrangements and possibly droplet rotation. This 
behavior can be contrasted with the complex relaxation 
of gels observed by Trappe and coworkers [25| . Here q- 



dependent and q-independent relaxations were observed 
simultaneously i.e. the intensity autocorrelation function 
exhibited two separate shoulders: one with a Tc that var- 
ied with q and the other with a different Tc which did 
not vary with q. Hence separate ballisticlike and large- 
scale rearrangements were taking place within these gel 
samples. In our study the creaming and ballisticlike rear- 
rangements are combined in a single relaxation process. 

The transition to the late stage is characterized by a 
marked change in the sample dynamics (possibly corre- 
sponding to a transition from the dispersed droplet region 
to a dense cream pJi]). At the simplest level the sample 
relaxation slows down; however there are four other fea- 
tures involved. First, the q-independent aspect of the 
dynamics comes to an end. This is evident in the de- 
cay rate versus wave vector graph meeting the y-axis at 
the origin (see inset to Fig. [3)Jb)) and implies that the 
emulsion is becoming compact. Second, for these droplet 
sizes and volume fractions we find that the variation of 
the decay rate with wave vector is independent of the 
orientation of the wave vector; such isotropy is common 
to high volume-fraction granular systems settling under 
gravity [3l|. (In studies of sedimentation velocities us- 
ing direct imaging of particles at lower volume fractions, 
$ < 0.2, show sustained anisotropy in the velocity and 
its variance [12] [l^-) Third, the exponent /3 which de- 
scribes the dynamic structure factor falls to one or be- 
low. Fourth, the dynamics appear to be characterized by 
a large variance which peaks at the characteristic time of 
the system relaxation. Hence the system combines ballis- 
ticlike dynamics with intermittency (which has been ob- 
served previously for various systems) with a stretched 
exponential intermediate scattering function (which is 
unprecedented). The growth in the variance of fluctu- 
ations is likely to be a consequence of the increasingly 
correlated behavior of the dense cream. 

Often, soft glassy systems have dynamics that slow 
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down with age. Studies have been carried out for dense 
colloids, shaving foam, surfactant 'onions' and many oth- 
ers. Cipelletti et al. [l^] and others [11] suggest that 
the characteristic dynamics are due to the build-up of 
internal stresses. The slowing down is then due to ei- 
ther a decrease in the rate of change of the stress source 
strength or a reduction in the number of active stress 
sources as the sample ages. In our case it is possible 
that the early regime is initiated either by formation or 
by the loading of the sample; this is strongly supported 
by the fact that it can (sometimes) be re-induced by re- 
moving the sample and replacing it again later. Hence, 
the two regimes we observe in the dynamics of particle- 
stabilized emulsions share much in common with the two 
regimes observed by Guo and coworkers in the dynamics 
of nanoemulsions 18]. Nanoemulsions are stabilized us- 
ing conventional surfactants and have subunits of much 
smaller size. It was observed using XPCS that during the 
period immediately following the loading of the sample 
the dynamics were fast; however, this behavior was ob- 
served to give way after 10-20 minutes to a regime char- 
acterized by slower dynamics. The initial fast dynamics 
were interpreted as being due to an initial population 
of stress dipoles created during the loading of the sam- 
ple; the drive behind the dynamics in the slower regime 
was not unambiguously identified. In our case the initial 
phase is likely to have dynamics that are driven both by 
the population of stress dipoles created during loading 
and also by the creaming under gravity. The drive due 
to gravity will not disappear in the second regime that 
we observe and this is likely to give rise to the late-stage 
dynamics. 

The linear g-dependence of 1 /tc such as that observed 
here in both regimes (Fig. [3]) has been attributed to 
stress dipoles within samples. It implies that the direc- 
tion of rearrangements within the sample is maintained 
across several rearrangements. Considering rearrange- 
ment events as randomly distributed in time Duri and 
Cipelletti [l^l have shown that if they are dipolar 
stresses then this necessarily leads to a /3-exponent of 
1.5 for low wave vector that approaches 1.0 for large 
wave vector. Our observations yield the same trend in 
(3 in the early stage; although, all of our measurements 
are carried out at large wave vector. For the late stage 
where intermittent events appear to dominate the stress 
source is still gravity pushing the particles upwards into 
each other. This no longer gives /?= 1.5. In all experi- 
mental observations to date the intermediate scattering 
function accompanying the ballisticlike dynamics could 
be described using a compressed exponential where the 
compression exponent was a little larger than 1. Re- 
cent experimental studies have shown that the dynamics 
are characterized by intermittent events [l^l [IB]. The 
results presented here differ in one key respect: the late- 
stage dynamic structure factor is stretched rather than 



compressed. Hence we have observed experimentally that 
ballisticlike dynamics together with large intermittent re- 
arrangements do not necessarily have a faster than expo- 
nential relaxation. Stretched exponential behavior does 
not emerge from the model of Bouchaud and Pitard [13] . 
An intermediate scattering function with stretched ex- 
ponential form is often observed for colloidal gel samples 
especially for measurements at low wave vector [34] . Here 
the dynamics are associated with structural heterogene- 
ity on the scale of clusters of colloids. In our experiments 
the dynamics have been characterized on the particle 
scale and so the same explanation is unlikely to apply. 
The cornpressed exponential relaxation of Bouchaud and 
Pitard [23] was associated with a long-range elastic de- 
formation field; the model predicts that the exponent will 
fall at larger wave vectors (and this has been observed). 
However an exponent of less than 1 does not sit easily 
with this interpretation. 



V. SUMMARY 

We have carried out an XPCS study of colloid- 
stabilized emulsions as they cream and form a novel soft 
solid. The data show two distinct regimes and a cross- 
over between them: the early time regime is characterized 
by a creaming of the droplets combined with local rear- 
rangements. As the samples age a peak in the variance 
of the correlation function is observed (most clearly for 
$1, =5%) with the peak located at the characteristic re- 
laxation time. This strongly suggests that the droplets 
undergo intermittent rearrangements during compactifi- 
cation. The increasing dominance of large fluctuations 
suggests that large regions of the sample are becoming 
strongly correlated. The form of the intermediate scat- 
tering function is unusual for systems exhibiting ballisti- 
clike dynamics with an intermittent character. Large, 
colloid-stabilized, droplets form intriguing systems for 
further study: the ability to select the size of the irre- 
versibly trapped particles provides a means to control 
the roughness of granular scale objects. 
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